Optimal Quantum Phase Estimation 
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By using a systematic optimization approach we determine quantum states of light with definite 
photon number leading to the best possible precision in optical two mode interferometry. Our treat- 
ment takes into account the experimentally relevant situation of photon losses. Our results thus 
reveal the benchmark for precision in optical interferometry. Although this boundary is generally 
worse than the Heisenberg limit, we show that the obtained precision beats the standard quantum 
limit thus leading to a significant improvement compared to classical interferometers. We further- 
more discuss alternative states and strategies to the optimized states which are easier to generate 
at the cost of only slightly lower precision. 

PACS numbers: 03.65.Ta, 06.20.Dk, 42.50.Lc, 42.50.St 



Interferometry is one of the most important measure- 
ment techniques in physics. Its numerous variations in- 
clude Ramsey spectroscopy in atomic physics, optical in- 
terferometry in gravitational wave detectors, laser gyro- 
scopes or optical imaging to name but a few. All these 
applications aim to estimate the quantity of interest, nor- 
mally a relative phase gathered by one "arm" of the inter- 
ferometer, with highest possible precision. In this letter 
we present fundamental limits to this precision in optical 
interferometry for light with definite photon number in 
the presence of losses. To find these limits it is neces- 
sary to consider the "cost" of the experiment, i.e. the 
required resources, and determine the precision of the 
estimated phase as a function of the cost. In optical in- 
terferometry the required resource is typically identified 
to be the number of photons N necessary to reach a de- 
sired precision. Classically the precision of the estimated 
phase scales then like 1/Vn, the so called standard quan- 
tum limit (SQL). Quantum interferometry on the other 
hand promises to beat this limit by employing highly non- 
classical entangled states to drastically improve the preci- 
sion to a scaling 1/A^ known as the Heisenberg limit [lil2]- 
The realization of interferometric measurements beyond 
the SQL is a very active field and recent years have seen 
tremendous progress [3|, Ij, S a ll| ■ A quantum enhance- 
ment in precision would allow for a significant reduction 
of the energy flux while keeping the same measurement 
precision. This is importantj^ for example, if the phase is 
induced by a fragile sample [3| . However, most of the the- 
oretical work done so far ignores the unavoidable presence 
of noise in the system. Existing treatments come to the 
conclusion that the benefit from highly entangled states 
deteriorates quickly even if only a small amount of noise 
is present in the system [8|, M, lid lul, ll4l ■ This is not re- 
ally a surprise since states of this type are typically very 
fragile: In optical interferometry, the well-studied NOON 
state promises to provide Heisenberg limited sensitivity 



in phase estimation [il3] , however, the loss of merely a 
single photon renders this state useless since it collapses 
into a product of two Fock states which cannot acquire 
any phase information. 

The Heisenberg limit is believed to be the ultimate 
precision in optical phase estimation, however, it is yet 
an unsolved problem if this limit can be reached in the 
presence of noise and, if not, then what is the ultimate 
precision? In this letter we answer this question for op- 
tical two-mode interferometry in the presence of photon 
losses, which is the limiting source of noise in such ex- 
periments. By using a systematic approach we determine 
optimal states with definite photon number leading to the 
highest possible precision. Although it turns out that the 
Heisenberg limit is unattainable we show that one can 
beat the SQL thus greatly improving precision beyond 
classical interferometry. Furthermore we introduce alter- 
natives to the optimal states, with simpler structure, at 
the cost of only slightly less precision. 




FIG. 1: General optical interferometric setup with two arms. 
Ghannel a acquires a phase ip relative to channel b. The two 
beam splitters in channel a and b symbolize photon losses. 

We consider a general interferometer with two arms as 
shown in Fig. [U in particular we do not put any restric- 
tions on the measurement scheme. Channel a is accu- 
mulating a phase tp relative to channel b and both arms, 
a and b, are subject to photon loss which can be seen 
as the effect of fictitious beam splitters inserted at ar- 
bitrary locations in both channels. We aim to estimate 



(/? with the highest possible precision quantified by using 
the uncertainty of the estimated phase ipest , 



Sip = ( (ipest\d{ipest)/dip\ -if 
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which, for an unbiased estimator, is simply the standard 
deviation. According to the general theory of quantum 
parameter estimation IJ, llSl . Ilq 5ip is bounded by the 
quantum Cramer-Rao bound 
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where v is the number of experimental runs and Fq is 
the quantum Fisher information. It was shown that this 
bound can always be reached asymptotically by maxi- 
mum likelihood estimation and a projective measurement 
in the eigenbasis of the 'symmetric logarithmic derivative 
operator' [ij]. Hence, inequality ^ defines the prin- 
cipally smallest possible uncertainty in phase estimation 
the determination of which is the primary scope of this 
paper. An explicit construction of the measurement op- 
erators will be given elsewhere [TTJ . 

Photon losses can be modelled by inserting fictitious 
beam splitters with transmissivities rja,b into both arms 
of the interferometer which couple the system to (un- 
correlated) environments. The noise operation and the 
phase accumulation commute, i.e. it is irrelevant if 
photons are lost before, during or after channel a ac- 
quires its relative phase with respect to b. If the noise 
operation is applied first, the state after tracing out 
the environmental degrees of freedom can be written as 
Z]m=o ^UaKkMPinKlj^K]^^ with Kraus operators 

= {^ — 'ila)^ila d'/vH, where a is the annihila- 
tion operator for mode a, and analogously for mode 6. 
This state acquires a phase through the transformation 
p{(p) — e~"^" "pe"^'^ °. This scenario is equivalent to 
a continuous photon loss model described by a master 
equation with loss rates | ln77a,b| per unit time. 

We consider the most general pure input states with 
definite photon number N, 
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= y^ak\k,N~k), 



(3) 



where \k,N—k) abbreviates the Fock state \k)a\N — k)i,. 
Special cases of ([3]) comprise in particular the highly en- 
tangled NOON state, (|7V0) + \0N))/V2 which, in the ab- 
sence of noise, leads to Heisenberg limited phase esti- 
mation but is strongly prone to decoherence otherwise. 
Equation ([3]) includes states which are "less" entangled 
but more robust representing a trade-off between preci- 
sion and robustness. Also, in the absence of additional 
reference beams, a superposition of states with different 
definite photon number would effectively become a mix- 
ture 18| of these states. Convexity of Fq 'l^ implies 



then that the analysis can be restricted to states with 
definite photon number if we use them successively [17 1 . 
In the case of no losses, the state of the system, 
\^{ip)) = e~"^° "|'0)m, remains pure and the quantum 
Fisher information reads 



Fq = 4[(V/MI^'M> - K^'MIVM)!'] = 4(A(ata))2, 

(4) 
where (A(a^^a))^ is the variance of a^a with respect to 
the state \^) in a nd the prime indicates a derivative with 
respect to ip |14l | . In the presence of noise the pure input 
state will deteriorate into a mixture p{ip). If the eigen- 
values and eigenstates of p{ip) are known th e q uantum 



Fisher information can be easily calculated [lj|. How- 
ever, very often the analytical diagonalization of p{ip) 
turns out not to be feasible. In this case, if the density 
operator is given in the form p{(p) = J2 Pjl^j iv)) {'^j {v) 1 1 
where the \ipj{p)) are not necessarily orthogonal, we can 
use the convexity of -Fq [ISl to obtain an upper bound 
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where the variance corresponds to |f/'j((/j)). The bound is 
reached if the spaces spanned by {|V'j(<p)), |V'j(<p))} and 
{|V'i(</5)), IV'U'/'))} are orthogonal for j 7^ «■ Particularly, 
we have Fq = Fq for the NOON state and, generally, if 
photon losses are only present in one channel, i.e., rjb — 1. 
The latter is relevant if the phase ip is induced by a sample 
in arm a which also causes the majority of photon losses. 
Applying Eq. ([5]) to the state ([3]), we get 
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For i]b = 1 we have 
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Obviously, the phases of the Uk are 
irrelevant. Furthermore, we proved analytically that -Fq 
and Fq are concave functions of the {xk} [iZl]- This 
simplifies the numerical maximization of Fq or Fq and 
more importantly, it implies that any maximum is global. 
Figure [H shows the results of such an optimization for 
Va = rjb = V = 0.9 and T]a = r] ^ 0.9, r/b = 1, i.e. 10% 
losses in both arms and one arm, respectively (blue, solid 
lines). In the following we concentrate on these two sce- 
narios. The quantity we analyze is S(pmin = I/^/Fq 

(or 1/ \ Fq) corresponding to the best measurement pre- 
cision for fixed v [see Eq. ([5])]. The lower and upper 
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FIG. 2: Phase estimation precision Stprnin for losses in both 
arms of the interferometer (a) versus photon number A*' (rja = 
rjb = r) = 0.9) and (c) versus transmissivity rja = rjb = rj 
{N — 20). The precision for losses in one arm, i.e. -qt — 1, 
is shown in (b) versus A^ {rja = rj — 0.9) and in (d) versus 
rja = rj {N = 20). Blue, solid line: Optimal state; Red, 
dashed hue: NOON state; Green, dashed-dotted line: NOON 
chopping strategy; Shaded area: Region between Heisenberg 
limit and classical limit (see text). 



boundaries of the shaded regions in Fig. [5] are the Heisen- 
berg limit, 1/iV, and a standard interferometric limit 
(SIL) |21| given by l/^/Wrj (losses in both arms) and 
(1 + ^)/2y/Nri (losses in one arm). Since the SIL is ob- 
tained by a classical reference experiment (it scales like 
the SQL), Sifmin falling into the shaded region implies an 
improvement over a classical interferometer. For rjt j^ 1 
we used the state which maximizes Fq to calculate Fq 
which differed by no more than 0.45%. Due to Eq. ^ the 
"true" maximum has to lie in between these quantities 
and its deviation from Fq can be neglected on the scale 
given by the difference of the SIL and the Heisenberg 
limit. As can be seen in Figs.[2la)(b) it is obviously not 
possible to reach the Heisenberg limit using input states 
with definite photon number. However, we gain a signif- 
icant improvement over the SIL of up to 60% (losses in 
both arms) and 73% (losses in one arm). The precision 
for the NOON state in the presence of losses [dashed lines 
in Figs.ma)(b)] is given by Sipmin — l/Nr]'^^^ for losses 
in both arms and ceases to be optimal for TV > 7 photons. 



For losses in one arm we have Sipmin = vl 
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which is generally worse than the optimal state: Here it is 
beneficial to use NOON states with unequal amplitudes of 
the two components. The best precision of such an "un- 
balanced" NOON state is given by (1 -I- rj-'^^'^)/2N which 
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FIG. 3: Coefficients Xk = \ak\ of the optimal state versus 
photon number A'^. (a) losses in both arms (rja = rjb = rj = 
0.9). (b) losses in one arm {rja — rj — 0.9, rjt = 1). 



coincides with the optimal state for N < 10. However, 
for larger photon numbers NOON states are not prefer- 
able; the precision gets even worse than the SIL. 

Figure [3] shows that the optimal state for losses in both 
arms has generally many non-zero components. Intu- 
itively this is consistent with the idea that the loss of 
a photon does not radically change the photon number 
distribution. The structure of the optimal state is sim- 
pler for losses in one arm. We therefore compare it to 
states with only two non-zero components. The best pre- 
cision obtained by these states differs by no more than 3% 
from the optimal case for the example shown in Fig.[2I^b). 
They have the form y^\m, N — m) + ^/l — p\N, 0) and are 
thus similar to the optimal state. This reflects the fact 
that it is both beneficial to have a large photon num- 
ber difference between arm a and b and have m > so 
that loss of a photon does not completely destroy the 
coherence. For equal losses in both arms, the best two- 
component state has approximately a symmetric form 
{\m,N — m) + \N — m,m))/^/2, but the corresponding 
precision deviates significantly from the optimal state [up 
to 13% for the example shown in Fig.I^^a)]. Here, states 
with more non-zero components are more useful, e.g. a 
twin Fock state [20| reducing the difference to 9%. Fig- 
ure [2jc) and (d) shows the best possible precision versus 
77 for iV = 20. For rj > 0.95 « e"^/^ the optimal state, 
the optimal two-component state and the (unbalanced) 
NOON state are identical. 

We can also use different strategies to operate the in- 
terferometer: Since our resources are given by the total 
number of photons, N, we can, instead of employing a 
single NOON state, use these photons to generate N/n 
nOOn states containing n < N photons each, i.e. we split 
up a "larger" NOON state into a number of "smaller" nOOn 
states which are sent successively through the interfer- 
ometer. Maximization over n (treated as a continuous 
parameter) leads to the precision 
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FIG. 4: Differential scaling of phase estimation precision 
Stpmin with number of photons A^ for different transmissiv- 
ities. (a) Optimal state for losses in both arms (rja = i]i, = i]) 
and (b) for losses in one arm {-qa = -q, rjb = 1). The color 
coding is the same in both plots. 



mentally generated. So even if the curves drops to 0.5 for 
N —f oo, interferometry using "smaller" quantum states 
has a significant advantage both in terms of scaling and 
absolute precision over classical interferometry. Partic- 
ularly for very small numbers of photons (unbalanced) 
NOON states are optimal. Above this threshold the use 
of the more complex, optimized states or one of our ex- 
perimentally more feasible alternatives is favourable. 
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where fj ^ fjo = 1, rjo — e and fj — -q, fjo = rjQ ^i 4.386 
for losses in both arms and one arm, respectively. In 
the latter case we use unbalanced nOOn states. Exam- 
ples for this "chopping" strategy are given by the green, 
dashed-dotted lines in Fig. [2l Note that in this case the 
total number of experimental runs is uN /n. The num- 
ber of photons "per run" is n — 1, n = ln?7o/| ln?7| and 
n = TV for the three cases in Eq. ([8]). For sufficiently 
small total photon numbers [cf. last line in Eq. ([8|)], it 
is not an advantage to chop the NOON state. If N is 
larger, the strategy does not improve the scaling with 
N compared to the SIL (or SQL). Nonetheless, it is an 
improvement over the SIL by a constant factor of approx- 
imately 2 (losses in both arms) or 2.5 (losses in one arm) 
in the example shown in Fig. [21 i.e. we need almost four 
(6.2) times less photons to reach the same measurement 
precision. 

The scaling of the precision of the optimal state with 
the number of photons turns out not to behave exactly 
like a power law. Therefore we define a differential scaling 
S{N) given by the local slope of 6ipmin{N) on a log- log 
scale obtained by a linear fit to the points corresponding 
to iV — 4, • • • , TV + 4. If Sipmin scales like a power law S 
would be the constant power, e.g. S = 0.5 for the SQL 
and S = I for the Heisenberg limit. Results are shown in 
Fig. m It is clearly visible that the scaling of the optimal 
state drops quickly, tending away from the Heisenberg 
limit towards the SQL. Even for rather high transmissiv- 
ity (green fines correspond to 95%) a scafing of ~ 0.61 
(losses in both arms) or '^ 0.68 (losses in one arm) for 
A^ = 70 is not exceeded. Moreover, the scaling gets worse 
for higher photon numbers. As yet it remains an unre- 
solved but fascinating question if the scaling of the op- 
timal states eventually tend to the SQL for all 77 < 1 
or asymptotically reaches a value which beats the SQL 
proving a true quantum scaling advantage in the presence 
of losses. Of course these are rather theoretical consider- 
ations: In practice there are restrictions to the size of the 
state (in terms of photon number) which can be experi- 
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